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Abstract 

We reformulate the calculation of the dark matter-nucleon scattering cross sec¬ 
tions based on the method of effective field theories. We assume that the scatterings 
are induced by the exchange of colored mediators, and construct the effective the¬ 
ories by integrating out the colored particles. All of the leading order matching 
conditions as well as the renormalization group equations are presented. We con¬ 
sider a Majorana fermion, and real scalar and vector bosons for the dark matter and 
show the results for each case. The treatment for the twist-2 operators is discussed 
in detail, and it is shown that the scale of evaluating their nucleon matrix elements 
does not have to be the hadronic scale. The effects of the QCD corrections are eval¬ 
uated on the assumption that the masses of the colored mediators are much heavier 
than the electroweak scale. Our formulation is systematic and model-independent, 
and thus suitable to be implemented in numerical packages, such as micrOMEGAs 
and DarkSUSY. 



1 Introduction 


Weakly Interacting Massive Particles (WIMPs) have been widely regarded as the most 
attractive candidate for dark matter (DM) in the Universe. They are weakly coupled to 
the Standard Model (SM) particles so that they are thermalized in the early Universe. 
Their relic abundance is determined by their annihilation cross sections, and it turns out 
that WIMPs may have a correct value of the cross sections to give the observed DM density 
Cl DM h 2 = 0.1196 ±0.0031 0 Moreover, such particles are often predicted in new physics 
beyond the SM. For instance, the lightest neutralino in the minimal supersymmetric 
Standard Model (MSSM) is a well-know candidate for WIMP DM. 

Since WIMPs are interacting with the ordinary matters, it is possible to use these 
interactions to detect them directly. Experiments based on such a strategy are called the 
DM direct detection experiments. These experiments search for the scattering signals of 
WIMPs kicking off target nuclei on the earth by detecting the recoil energy transferred to 
the nuclei. At present, the LUX experiment has the best sensitivity, and provides a limit 
on the spin-independent WIMP-nucleon scattering cross section as Ogi < 7.6 x 10” 46 cm 2 
at a mass of 33 GeV |2|. Further, there exist several proposals with ton-scale detectors, 
which will offer extremely improved sensitivities. 

To study the properties of WIMPs based on the direct detection experiments, it is 
necessary to evaluate the WIMP-nucleus scattering cross sections accurately. The in¬ 
teractions of a WIMP with a nucleon, as well as a nucleus, are generated through the 
couplings of the WIMP with quarks and gluons. These couplings are described in terms 
of the parton-level interactions, and contribute to the interactions with a nucleon through 
non-perturbative QCD effects. An appropriate way to compute the contribution is to take 
an effective theoretical approach. Here, the parton-level interactions are expressed by the 
higher-dimensional operators, and their contribution to the WIMP-nucleon couplings is 
computed by means of their nucleon matrix elements. See Refs. [3}j4] for the treatments. 

Usually, the parton-level interactions are mediated by heavy particles. Since the 
WIMPs are singlet with respect to the SU(3)c®U(1)em symmetry, the particles medi¬ 
ating the WIMP-quark interactions should be also electrically neutral and color singlet 
when they are exchanged in the f-channel, while they should be charged and colored when 
they are exchanged in the s- or ?j- channel. The extent of the significance of these contri¬ 
butions highly depends on models. For example, in the case of the neutralino DM in the 
MSSM, the Z boson and the Higgs boson mediating processes are classified into the for¬ 
mer type, while the squark exchange is the latter one. When squarks are extremely heavy, 
only the former contribution is sizable. In the limit of pure gaugino or higgsino case, on 
the other hand, the former contribution vanishes and thus the latter may be dominant. 
Therefore, it is desirable to construct a formalism to evaluate all of the contributions 
precisely enough on an equal footing. 

Recently, the LHC experiments give stringent limits on the masses of new colored 
particles. The results may suggest that the colored mediators which induce the couplings 
of WIMPs with quarks and gluon, if they exist, should have masses much heavier than ~ 
100 GeV. In this case, the contribution of such a particle to the WIMP-nucleon interaction 
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receives sizable QCD corrections because of the mass hierarchy and the large value of the 
strong coupling constant. This motivates us to reformulate the calculation in the following 
way; first, we obtain the effective theory which consists of the higher dimensional operators 
of DM and quarks/gluons at the energy scale of the mediator mass. At this point, we need 
to match the effective theory to the full theory so that only the short-distance contribution 
is to be included in the Wilson coefficients of the higher-dimensional operators. Then, 
we evolve the operators by using the renormalization group equations (RGEs) down to 
the scale at which the nucleon matrix elements of the operators are evaluated. It is the 
prescription that we will discuss in this paper. We will formulate a complete framework 
to carry out the calculation to the leading order in the strong coupling constant. In this 
formulation, all of the model dependence is included into the Wilson coefficients of the 
effective operators, and the rest part of the steps is independent of models. Therefore, 
the method is quite suitable for generic computational codes such as raicrOMEGAs |h] or 
DarkSUSY |6], 

In addition, we will discuss in detail the treatment for the twist-2 operators and show 
that we do not have to evolve the operators down to the hadronic scale (~ 1 GeV); we may 
evaluate their nucleon matrix elements around the electroweak scale. This point is often 
misunderstood, and we believe that our analysis clarifies the confusion. Furthermore, we 
will show the significance of the renormalization effects on the operators when the masses 
of the colored mediators are much heavier than the electroweak scale. 

This paper is organized as follows. In Sec. [2j we discuss our formulation in the case 
where DM is a Majorana fermion. We present the effective operators for the Majorana 
fermion and evaluate the matching conditions on their Wilson coefficients. The RGEs of 
the operators are listed there as well. The results for the real scalar and vector boson 
DM cases are also given in Sec. [3] and Sec. [4j respectively. Then, in Sec. [5j we study the 
renormalization effects by using a particular model. Also, we will discuss the treatment 
of the twist-2 operators in the section. Sec. [6] is devoted to conclusion. In Appendix, we 
present the formulae for the one-loop contribution to the scalar-type gluon operators in 
the case of scalar DM, which as far as we know have not been given in the literature so 
far. The result is useful when the colored mediators have masses around the electroweak 
scale. 


2 Formalism: Majorana fermion DM 

In this section, we give a formalism to evaluate the scattering cross sections of a WIMP 
with a nucleon. The procedure described here consists of the following steps. First, we 
construct the effective theory for the WIMP, quarks, and gluons, by integrating out the 
mediator particles. The effective interactions obtained here are expressed in terms of the 
higher-dimensional operators. Then, we evolve the Wilson coefficients of the effective 
operators according to the RGEs down to the scale at which the nucleon matrix elements 
of the operators are evaluated. Finally, we express the effective coupling of the WIMP 
with a nucleon in terms of the Wilson coefficients and the nucleon matrix elements. The 
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scattering cross sections are readily obtained from the effective coupling. We will evaluate 
them to the leading order in the strong coupling constant throughout this work. 

In this section, we assume that the WIMP DM to be a Majorana fermion. The 
real scalar boson DM and the real vector boson DM cases are discussed in Sec. |3] and 
Sec. [4J, respectively. It should be noted that Dirac fermion and complex scalar DM 
candidates are severely constrained by the DM direct detection experiments because they 
in general couple to the vector current of quark fields. Furthermore, this vector interaction 
is not renormalized, and thus the conventional way of calculation is sufficient for this 
contribution |7|. For these reasons, we do not consider these cases in this paper. 


2.1 Effective Lagrangian 

To begin with, let us write down the effective interactions of a Majorana fermion, which is 
assumed to be a WIMP, with quarks and gluon. The interactions are expressed in terms 
of the following higher-dimensional operators |3|: 


= E + E E C T,° r T, + E C AvO\v . (1) 

p=q,g i=\ .2 P=q,g q 

with 

Of, = , 

Of = -F X°G^G a ^ , 

71 

o? 1 = . 

°r, = , 

°av = • ( 2 ) 

Here, we only keep the operators that remain sizable in the non-relativistic limit. In 
addition, we have used the classical equations of motion and the integration by parts to 
drop the redundant operators |8,9j|. y°, q, and G^ )V denote the Majorana fermion, quarks 
(q = u,d, s,c,b,t), and the field strength tensor of gluon field, respectively; m q are the 
masses of quarks; M is the mass of the WIMP; a s = g 1 2 s /{ 47t) is the strong coupling 
constant, 0^ v and O 9 ^ are the twist-2 operators of quarks and gluon, respectively, which 
are defined bjj^] 


= \qi+ D vi„ - jvE i . 

oy = Gfai P - , (3) 


1 Notice that we have changed the definition of 

follow the convention used in Ref. 


by a factor of —1 from those in Refs. 3 4 . We 
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Tabic 1: Mass fractions. These values are based on the lattice QCD simulations 


13 


14 


Proton 

Neutron 

/£? 0.019(5) 
0.027(6) 
0.009(22) 

.ft ] 0.013(3) 
fS 0.040(9) 
0.009(22) 


with the covariant derivatives. The effective operators are defined at the mass scale 
of mediators, which is assumed to be well above the mass of top quark. Generalization to 
other cases is straightforward; for instance, if such heavy particles have masses similar to 
or lighter than the top mass, one should integrate top quark as well so that the effective 
theoretical approach is appropriate 0 

Note that we include a s /ir to the definition of the gluon scalar-type operator 0 9 s . We 
discuss the meaning in the next subsection. 


2.2 Nucleon matrix elements 


As discussed in Introduction, we need the nucleon matrix elements of the effective opera¬ 
tors to evaluate the WIMP-nucleon effective coupling. These operators are classified into 
three types in terms of the Lorentz transformation properties of the quark bilinear parts 
in the operators; the scalar-type operators (Og, 0 9 s ), the axial-vector operator (0\ v ), and 
the twist-2-type operators (PfQ. Since these operators do not mix with each other 

under the renormalization group (RG) flow, we are allowed to consider them separately. 

As for the scalar-type quark operators O q s , we use the results from the lattice QCD 
simulations. The expectation values of the scalar bilinear operators of light quarks be¬ 
tween the nucleon states at rest, | N) (N = p, n), are parametrized as 


/f } - 
J 1( 1 


= (N\m q qq\N)/m N , (4) 

which are called the mass fractions. These values are shown in Table [lj Here, is the 


nucleon mass. They are taken from Ref. 
by using the results from Refs. 


13,14 


12 , in which the mass fractions are computed 


The nucleon matrix element of 0 9 s is, on the other hand, evaluated with the trace 


anomaly of the energy-momentum tensor 15 . For Nf = 3 quark flavors, the trace of the 


energy-momentum tensor in QCD is given as 

= _^ G a G a^ 
m 8vr ^ 


Y m i qq > 

q=u,d,s 


(5) 


2 In Ref. 
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such a situation is discussed where the exchanged particle has a similar mass to the 


6-quark mass. In this case, of course, 6-quark (also top quark) should be simultaneously integrated out 
when the effective theory is formulated. 
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Tabic 2: Second moments of the PDFs of proton evaluated at p — mz • We use the CJ12 
next-to-leading order PDFs given by the CTEQ-Jefferson Lab collaboration 


16 


9(2) 

0.464(2) 



«(2) 

0.223(3) 

8(2) 

0.036(2) 

d(2) 

0.118(3) 

8(2) 

0.037(3) 

s(2) 

0.0258(4) 

8(2) 

0.0258(4) 

c(2) 

0.0187(2) 

8(2) 

0.0187(2) 

b(2) 

0.0117(1) 

6(2) 

0.0117(1) 


up to the leading order in a s . The relation beyond the leading order in a s is also readily 
obtained from the trace-anomaly formula. By evaluating the operator (J5]) in the nucleon 
states |IV), from {N | | N) = we then obtain 


(N^G^G^N) = -?.»„/« , ( 6 ) 

with fj?) = 1 — J2 q =uds Notice that the r.h.s. of Eq. (j6j) is the order of the typical 
hadronic scale, O(mjv). That is, although we include a factor of a s /Ti in the definition 
of 0 9 s , its nucleon matrix element is not suppressed by a s /7i. This is the reason why we 
have defined 0 9 s to contain Oi s ji\. 

Next, we discuss the nucleon matrix elements of the twist-2 operators. They are given 
by the second moments of the parton distribution functions (PDFs): 


with 


(N{p)\O q \N(p)) = -^-(p„Pv-\rn 2 N g fl J){q{2-,p) + q(2]p)) , 

Til n \ / 

( n (p)\°U n (p^ = - \ m N9n^)g& p) ■ 


q{ 2 ;f) 

g( 2;/x) 



x q(x, p) 



x q(x, l-i) 



x g(x,p) 


J 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Here q(x,p), q(x,p) and g(x,p) are the PDFs of quarks, antiquarks and gluon at the 
factorization scale p, respectively. These values are well measured at various energy 
scales, contrary to the case of the scalar matrix elements. In Table [2j for example, we 
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Figure 1: Tree-level matching condition for Majorana fermion-quark effective interactions. 
Gray square represents the vertex for the quark effective operators. 


present the second moments for proton at the scale of /i = mz with mz the Z boson 
mass. Here, we use the CJ12 next-to-leading order PDFs given by the CTEQ-Jefferson 
Lab collaboration [16]. Those for neutron are given with the exchange of up and down 
quarks. As can be seen, the second moment for gluon g( 2; /i) is of the same order of 
magnitude as those for quarks. As a result, the nucleon matrix element of the gluon 
twist-2 tensor in Eq. ([8]) is This justifies the definition of where we have not 

included a factor of a s /tt in this case. Our definition for the gluonic operators (0 9 s , and 
O a T ) clarifies the order counting with respect to a s /n. 

Finally, the nucleon matrix elements of the axial vector-type operators are given by 


(^kf7/FM|A0 = 2s^Aq N . 


( 12 ) 


with A(pv called the spin fractions and s^ being the spin of the nucleon. The values of 


the spin fractions are taken from Ref. 17 ; A u p = 0.77, A d p = —0.49, and A s p = —0.15 


for proton. Those of neutron are to be obtained by exchanging the values of up and down 
quarks. 


2.3 Wilson coefficients 

Next, we evaluate the Wilson coefficients of the effective operators by integrating out 
heavy mediator particles. Here, we consider a generic situation in which the interaction 
Lagrangian of the Majorana fermion with quarks is given by 

Ant = q{a q + 5 g 7 5 )x 0 9 + h-c. , (13) 

where q denotes a heavy, colored scalar particle with its mass represented by M q . 

The interaction gives rise to the coupling of the Majorana fermion with quarks via 
the tree-level exchange of the colored mediator. By evaluating the diagram, we readily 
obtain the Wilson coefficients of the effective operators containing quarks. The matching 
procedure is illustrated in Fig. [l] Here, the gray square represents the vertex for the 
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Figure 2: One-loop matching condition for Majorana fermion-gluon effective interactions. 
Black square represents the vertex for the gluon effective operators. 


WIMP-quark effective operators. As a result, we have 




Cfi (/4f) 
C'To^f) 
Cav^f) 


a 2 + b 2 q M a 2 q -b 2 q 1 

8 (Mf - M 2 ) 2 ~~ 4 m q Mf — M 2 ’ 

a 2 + b 2 M 

2 (Mf - M 2 ) 2 ’ 

0 , 

i 

4 Mf - M 2 ' 


(14) 


Here, fj,p denotes the factorization scale, which is taken to be the mass scale of exchanged 
scalar particles. We have performed the expansion in terms of the quark momenta; he., 
this calculation is valid when (p ■ q)/(M~ — M 2 ) <C 1 with p M and q y, being the momenta 


of DM and quarks, respectively (See also footnote | 

Next, we derive the matching condition for the WIMP-gluon effective operators. At 
this point, we need to include only the short-distance contribution to the Wilson coeffi¬ 
cients [4,19. This is achieved by the matching procedure shown in Fig. [2j This reads 


CIM = - 

C^_(pf) — 0 • 


M(a 2 + b 2 ) 

^ 96Mf (Mf - M 2 ) 


(15) 


3 On the other hand, it turns out that the WIMP-nucleon scattering cross sections are considerably 


enhanced when M q — M < 100 GeV 18 . 
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Here, the summation is taken over all quark flavors. Note that the Wilson coefficient of 
the gluon scalar operator is generated at the leading order in a s , he., 0[a °), while those 
of the gluon twist-2 operators vanish at this order |4|. They are induced at 0(a s ). 


2.4 Renormalization group equations 

The effective operators obtained above are evolved by means of RGEs. In this section, we 
list the RGEs for the operators which we use in the following analysis. In this paper, we 
only use the one-loop RGEs since our main concern is to formulate the procedure for the 
leading-order calculation in a s . For the sake of convenience, however, we also mention 
some results that may be used for the higher-order calculation^] 

The one-loop beta function of the strong gauge coupling constant is given by 


da s 



f3(a s ) 


2vr 



N c + ~ N f 


(16) 


where N c = 3 is the number of colors and Nf denotes the number of quark flavors in an 
effective theory. Higher-order contribution is found in Refs. 


23 


Now we give the RGEs for the Wilson coefficients of the above operators. First, we 
consider the RGEs for the scalar-type operators (Of, Of). To that end, notice that the 
quark mass operator is RG invariant in a mass-independent renormalization scheme like 
the MS scheme, i.e., 

fi—m g qq = 0 . (17) 

d/i 

Then, by differentiating the trace anomaly formula ([5]), we also find 




(18) 


as Eq. ([5]) is an operator equation and thus scale-invariant. Accordingly, the scalar-type 
operators are RG invariant at 0(a s ). This is another reason why we include a factor of a s 
in the definition of Of. Beyond the leading order, the scalar-type gluon operator runs and 
mixes with the scalar-type quark operators during the RG flow. The RGEs for the case 
are obtained again by using the trace-anomaly formula with the use of the higher-order 
beta function of the gauge coupling constant 22,23 and anomalous dimensions for quark 
masses 24,25]. 

Next, we consider the RGEs for the twist-2 operators ( O f., Of(). The one-loop anoma¬ 
lous dimension matrix of the operators is evaluated as 1261 


(19) 


1 ! t ■ ( r 1 — ' tVrJ Tr , 


4 See also Refs. 20 21 for relevant discussion. 














with r T a (Nf + 1) x (Nf + 1) matrix: 


/¥ c F 


a. 


Tt ~^ 


0 

3 


4 \ 

3 






o 


lb 

~3 C ' f 


4 

3 F 3 

16 ^ *N f J 


( 20 ) 


where Cp = 4/3 is the quadratic Casimir invariant. Higher order RGEs are found in 


Ref. 27 


Finally, the RGE for the axial-vector interaction is readily obtained since at the leading 
order the axial-vector current is conserved: 


d „„ 

^ Cav = 0 ■ 


( 21 ) 


For higher-order corrections, see Ref. 28 


2.5 Quark threshold matching 

During the RG flow, another matching procedure is required when one goes across a 
quark threshold. For instance, around the t-quark mass threshold /r 4 ~ m t , the Wilson 
coefficients are matched as 


QH/h) IiVf=5 — Cs(/h )\ Nf =6 , 
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1 + — a s (fi t ) 


C 9 s (nt)\N f =5 - ’—^2 

C ! T 1 (^\ N f =5 = C^.(nt)\N f = 6 , 
C ! T l (^t)\Nf=5 = C^.(Ht)\N f =6 , 
^Avi^t) |n /=5 = C'AvifJ't) |iV /=6 ) 


Qs(Ah)|iV/=6 + Cl(/ i 0lA r /=6 j 


( 22 ) 


with q = u,d, s, c, b. Notice that although we consider the leading-order calculation, we 
include the next-to-leading order constant contribution to C 9 S here, since the effect is 
known to be large 29 . Similar matching should be carried out at the b- and c-quark 


threshold scales, if the coefficients are evolved down below them. 

In addition, at the threshold scales, the higher-dimensional operators suppressed by a 
power of the corresponding quark masses may also be generated. For example, at a heavy 
quark mass threshold rriQ, the ijiqQQ operator gives rise to not only the scalar-type gluon 
operator —a s (mQ)G^G j4 ^/(127r), but also the following dimension-six operators 


30 


31 


OIs^FIq) , A W nP( ^A -v 9 sOis(mQ ) A rlB n pr ,C 

60tt m 2 n D Gv ^ D G pf> 79rwrr.2 JabcG^G G vp , 


720tT771q 


(23) 
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with f abc th e structure constant of SU(3)c- Among them, those induced at the charm 
quark threshold may yield a significant effect. The naive dimensional analysis tells us 


that the higher dimensional operator (23) might give a correction up to AQ CD /m 2 ~ 10% 
to the leading term, though the correction may be parametrically suppressed to a few % 
by the prefactors of those operators. 

At present, there is no way to estimate the contribution of the higher-dimensional 
operators more accurately. Thus, it should be considered as a theoretical uncertainty of 
the computation. One way to reduce the uncertainty is to use the nucleon matrix elements 
evaluated above the charm threshold. As will be discussed in Sec. [5j for the twist-2 
operators, it is possible to use the PDFs evaluated above the charm/bottom threshold. 
For the scalar-type operators, on the other hand, the present lattice simulations are not 
able to precisely evaluate the charm-quark matrix element 32 . Future simulations may 


compute it with sufficient accuracy and help to reduce the theoretical uncertainty. 


2.6 Scattering cross sections 

Finally, we obtain the effective coupling of the Majorana fermion with a nucleon. For the 
spin-independent coupling, we have 


41° = InX^NN , 


(24) 


where N denotes the nucleon field and 


3 

4 


f N /m N = Y C'K^had )f? q ] - gC'K^had)/^ 
q=u,d,s 
Nf 

EE CriMtoC 2 ;/*) + ?( 2 ;/^)] - Y c tM) 9 &iA , (25) 

q 4 = 1,2 4 = 1,2 

where /ihad is the hadron scale usually taken to be around 1 GeV. For the contribution 
of the twist-2 operators, the summation runs over the number of the active quark flavors 
in the effective theory at the energy scale of /i where their nucleon matrix elements are 
evaluated. As will be shown in Sec. [5j the scale fi does not need to be taken at the 
hadronic scale; it is allowed to be set around the electroweak scale as long as the PDFs 
of the scale are known. 

The spin-dependent effective coupling is, on the other hand, given by 

£se? = a ivX°7 M 75X°A", (26) 

with 

a,N = E C q AV Ag N . (27) 

q=u,d,s 

By using the effective couplings, we finally obtain the scattering cross section of the 
Majorana fermion with a target nucleus as follows: 

2r _ , 

(28) 


4 

a = — 


MMt 


n \ M + Mt 


I Upfp T IT'nfnl'* T 4 T |flp(Sp) T On(s n )| 


J 
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where Mt is the mass of the target nucleus; n p and n n are the numbers of protons and 
neutrons in the nucleus, respectively; J is the total spin of the nucleus, and (sn) is the 
expectation value of the spin of a nucleon in the target. Here, we calculate WIMP- 
nucleus scattering cross sections in the limit of zero momentum transfer, for which each 
WIMP-nucleon scattering amplitude adds up coherently [7]. 

3 Formalism: real scalar boson DM 

Next we briefly show the results for the case of real scalar boson DM. We may use a 
similar procedure to that given in the previous section to formulate effective theories for 
the WIMP. 


3.1 Effective Lagrangian 

The effective interactions of the real scalar </> with quarks and gluon are expressed by 


with 


c, s = V dO’ s + V Cf, , 

p=q, 9 p=q,g 


(29) 


Of = tt> 2 m,qq , 

Of = -<fG A ^G% . 

7r 

o% = T . (30) 

Note that there is no spin-dependent interactions in the case of scalar boson DM. 


3.2 Wilson coefficients 


We next discuss the matching condition for the Wilson coefficients of the above operators 
in a theory where the interactions of the scalar boson with quarks are given by 


£ = i>,(a q + b^qtf) + h.c. , 


(31) 


where i/j q denotes a colored fermion with a mass of M. Then, with the tree-level matching 
procedure for the WIMP-quark interactions illustrated in Fig. [3j we obtain 


Gs(^f) 


Gt 2 (^f) 


< + K 2M l- M ' 2 Mi,, 

2 (Ml - M 2 ) 2 + m. Ml - M 2 ’ 
2(a 2 + b 2 )M 2 
(Ml - M 2 ) 2 ■ 


(32) 

(33) 
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Figure 3: Tree-level matching condition for scalar boson-quark effective interactions. Gray 
square represents the vertex for the quark effective operators. 


Again, the calculation is valid only when the mass difference between the heavy mediator 
particle and the real scalar boson is much larger than the energy of external quarks. The 
matching condition for the gluon operators is, on the other hand, obtained through the 
procedure shown in Fig. [4j We have 


c sM E 12 (M^ - AP) ’ 

(34) 

C'f- 2 (/Jf) = 0 . 

(35) 


Further, in Appendix, we give a result for the loop-computation of the one-loop diagrams 
in the left-hand side in Fig. [4| since as far as we know there has been no such a calculation 
in the literature. The result is useful for the cases where the masses of the colored particles 
are not so heavy compared with some of the quark masses. 

3.3 Scattering cross sections 

We now ready to evaluate the scattering cross section of the real scalar boson with a 
target nucleus. The spin-independent coupling of the real scalar boson with a nucleon 
defined by 

4V = fNcfNN , (36) 

is evaluated as 

h/m K = Y. CsMfrP ~ ICsMfra 

q=u,d,s 
A If 

+ + ?( 2 ;/^)]- • (37) 

q 
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Figure 4: One-loop matching condition for scalar boson-gluon effective interactions. Black 
square represents the vertex for the gluon effective operators. 


In the scalar boson case, there is no spin-dependent coupling with a nucleon. By using 
the effective coupling, we calculate the scattering cross section of the real scalar boson 
with a target nucleus as follows: 


1 ( Mx 

7T \M + Mt 


Upfp Tlnfn 


(38) 


4 Formalism: real vector boson DM 


Finally, we consider real vector boson DM. For previous calculation, see Ref. 
references therein. 


33 and 


4.1 Effective Lagrangian 

The effective interactions of the real vector boson B^ with quarks and gluon are written 
as 

= £ CIO”, + £ c t 2 Ot 2 + Cav®av > ( 39 ) 

p=g,g p=q,9 <? 
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-3/1 B M 



3^ Bn 



Figure 5: Tree-level matching condition for vector boson-quark effective interactions. 
Gray square represents the vertex for the quark effective operators. 


with 

O q s = B^Bnm q qq , 

O a s = —B p B p G A ^G A v , 

o q T , = -B BndHa .- Bp0l , 

oV = , (40) 

where is the totally antisymmetric tensor with e 0123 = +1. Here, the vector boson 

field is supposed to satisfy the on-shell condition (□ + M 2 )H m = 0 and dpB^ = 0. 


4.2 Wilson coefficients 


Let us evaluate the Wilson coefficients of the above operators in the presence of a fermionic 
colored particle Tp q coupling to the WIMP and quarks through the interactions: 

jC = ^ q {a q Y + b q Ylb)qBn + h.c. . (41) 

In this case, Fig. [5] yields the matching condition for the WIMP-quark effective couplings 
as 


CUpf) 


_(„ 2 , .2,_ *, _ 

y " “ 2(M 2 q - M 2 ) 2 

2 K + gyg 

(Ml - M 2 ) 2 ’ 

*(■»? + 

VJ.-M 2 ' 


m, - M 2 ’ 


(42) 

(43) 

(44) 
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Figure 6: One-loop matching condition for vector boson-gluon effective interactions. Black 
square represents the vertex for the gluon effective operators. 


As for the gluon contribution, Fig. [6] reads 


Cs{MF) ^ 12 (Ml - M 2 ) ’ 

(45) 

C't 2 (l- 1 ?) = 0 • 

(46) 


4.3 Scattering cross sections 

By using the results obtained above, we now evaluate the spin-independent WIMP-nucleon 
coupling. With the definition 

4f = InB^NN , (47) 

we have 

f N /m N = ^ Of (/ihad)/^ - 

q=u,d,s 

A T f 

+ + 2(2;/^)] - \c 9 T2 (n)g{2-v) . (48) 

q 

On the other hand, the spin-dependent effective coupling is given by 

4d = B>N-i° ls N . (49) 
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DM mass: 200 GeV, scalar mass: 700 GeV 
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Figure 7: Comparison of twist-2 contributions to WIMP-proton effective coupling cal¬ 
culated with PDFs obtained at /i = 2 GeV and mz- Here, we assume the WIMP is 
a Majorana fermion coupled with right-handed t- and 6-quarks (see text), and we take 
M = 200 GeV and M q = 700 GeV. Red (light-pink) bar denotes the uncertainty coming 
from the PDF input (perturbation in a s ). 


where 

a N = ^2 C q AV Aq N . (50) 

q=u,d,s 

With these effective couplings, we eventually get the scattering cross section of the real 
vector boson with a target nucleus as 


1 

cr = — 




7t V M + Mt 


I Tlpfp V ^lnfn\ "F g j |Op(Sp) T On(^n) | 


(51) 


5 Analysis 

Now we apply our formulation to a concrete DM model, and discuss the renormalization 
effects on the calculation. Here, we consider a Majorana fermion interacting with only 
the third generation right-handed quarks with a unit coupling constant; i.e., a q = b q = 0 


for q = u,d,s,c and a q = b q = 1/2 for q = b,t in Eq. (13). Since the new colored 


scalar particles introduced to the scenario only couple to the third generation quarks, 
the LHC constraints on them are less severe compared with those interacting with the 
first two generation quarks. This model is a simplified model for a system composed of 
a neutralino DM and a pair of right-handed stop and sbottom in the MSSM, though the 
couplings between them are different from above assumption. 
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By using the model, we first discuss the scale at which we evaluate the nucleon matrix 
elements of the twist-2 operators. As mentioned in Sec. 2.2 the twist-2 operators are not 
mixed with the scalar-type and axial-vector operators so that we may take a scale for the 
matrix elements of the twist-2 operators which is different from that of the other ones, and 
the matrix elements are obtained in a wide range of energy scales. Thus, it is important 
to determine which scale is appropriate for the calculation of the twist-2 contribution. In 
Fig.0 we compare the twist-2 contributions to the WIMP-proton effective coupling, 

N f 

f p (twist-2) /m p = ^ EE ^(fM 2 ;^) + <?( 2;/i)] - > ( 52 ) 

q i=l,2 £=1,2 

evaluated with the PDFs obtained at // = 2 GeV and mz- Here, we take M = 200 GeV 
and Mq = 700 GeV. The red (light-pink) bar denotes the uncertainty coming from the 
PDF input (perturbation in a s ). For the estimation of the uncertainty from the PDF 
error, we follow the method described in Ref. 16 with the y 2 tolerance T taken to be 
T = 10. The uncertainty caused by the neglect of the higher-order contribution in a s is 
evaluated by varying the input and quark-mass threshold scales by a factor of two, i.e., 
Mq/2 < up < 2 Mq, m t /2 < /q < 2 m t , and so on. It turns out that both calculations 
predict similar values for the twist-2 contribution, though the theoretical error in the case 
of /i = 2 GeV is a little bit larger than that with /i = mz- If one sets /i = 1 GeV, 
we expect that the error becomes much larger due to the charm-quark threshold effects 
since in the low-energy region the strong coupling constant rapidly grows up. In addition, 
the higher-dimensional operators suppressed by a power of the quark masses may give 
significant contribution if the scale /i is taken to be at a low-energy scale as discussed 
in Sec. 2.5[ which also contribute to the theoretical uncertainty. For these reasons, we 
conclude that it is appropriate to set the PDF scale // in a high-scale region, not the 
hadronic scale. In the following calculation, we take the scale to be fi = mz- 

Next, we show the renormalization effects on the WIMP-nucleon scattering. As dis¬ 
cussed in Sec. 2.4, the twist-2 operators receive the renormalization effects. The effects 
are expected to be significant when the input scale \xf, be., the typical mass scale of 
colored mediators, is much higher than the electroweak scale. 

we show each contribution to the WIMP-proton effective coupling f p 


In Fig. 8(a^ 


as functions of the mediator mass Mq. We set the WIMP mass to be M = 200 GeV. 
The upper red (lower blue) line shows the contribution of the scalar-type (twist-2-type) 
operators. For the twist-2 contribution, we show both calculations with and without 
the renormalization effects in solid and dashed lines, respectively. By using the effective 
coupling, we then evaluate the WIMP-proton scattering cross section a p . We plot it as 
a function of Mq in Fig. 8(b) Here again the WIMP mass is set to be M = 200 GeV. 


The solid and dashed lines show the results with and without the renormalization effects, 
respectively. We find that the renormalization effects change the resultant value for the 
twist-2 contribution by more than 50% when Mq > 500 GeV. In this case, the scattering 
cross sections are modified by more than 20%. The results indicate that it is important 
to include the RGE effects, especially when the colored mediators are much heavier than 
the electroweak scale. 
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M=200 [GeV] 


M=200 [GeV] 




(a) Effective coupling f p (b) Scattering cross section cr p 

Figure 8: (a) Each contribution to the WIMP-proton effective coupling f p as functions 
of the mediator mass Mg. DM model adopted here is the same as Fig. [7J Upper red 
(lower blue) line shows the contribution of the scalar-type (twist-2-type) operators. For 
the twist-2 contribution, solid and dashed lines show the results with and without the 
renormalization effects, respectively, (b) WIMP-proton scattering cross section a p as a 
function of Mg. Solid and dashed lines show the results with and without the renormal¬ 
ization effects, respectively. In both plots, WIMP mass is set to be M = 200 GeV. 


6 Conclusion and discussion 


So far we have discussed a way of evaluating the WIMP-nucleon scattering cross section 
at the leading order in a s based on the effective theoretical approach. We have considered 
a Majorana fermion, real scalar and vector bosons, and presented formulation for each 
case. Further, using a particular example with a Majorana fermion, we have shown 
that the renormalization effects may change the twist-2 contribution by more than 50% 
when the colored mediators are much heavier than the clectroweak scale, which results in 
modification to the WIMP-nucleon scattering cross section by 0( 10)%. 

As shown in Fig. [Tj, the calculation of the twist-2 contribution suffers from (9(10)% 
uncertainty due to the perturbation in a s . It is possible to reduce the uncertainty by 
going beyond the leading-order calculation. In fact, we have already had the higher-order 
results for the RGEs and the matching conditions at each quark threshold, as commented 
in Sec. 2.4 To complete the next-to-leading order computation, however, we further need 


the higher-order matching conditions between the full and effective theories at the input 
scale. We defer the calculation as future work. In addition, we expect that future lattice 
QCD simulations will much improve the determination of the quark content in nucleon. 
These two developments will enable us to evaluate the WIMP-nucleon scattering cross 
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Figure 9: One-loop contribution to the WIMP-gluon coupling. 


sections with great accuracy. 

Finally, we would like to emphasize that the prescription for the computation of the 
WIMP-nucleon scattering cross section discussed in this paper is quite systematic, and 
the formulation itself is almost model-independent. The model-dependence is included in 
the Wilson coefficients at the factorization scale, and the subsequent procedure is similar 
in every case. Therefore, this method is suitable to be used in general computational 
codes for the direct detection rate of DM, such as micrOMEGAs [5] and DarkSUSY |6|. 
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Appendix: gluon-loop contribution for real scalar bo¬ 
son DM 


In this Appendix, we give a result for the calculation of one-loop gluon contribution in 
the case of real scalar boson DM, with quark masses kept non-vanishing. Similar results 
have been already obtained for the cases of Majorana fermion and real vector boson in 
Refs. |4| and |33 


respectively. 


We assume that the interactions of the real scalar boson with quarks and the cor¬ 
responding colored heavy fermions are described in terms of the Lagrangian given in 
(31). The diagrams we consider here are illustrated in Fig. [9} By evaluating the 


Eq 
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diagrams, we compute the contribution of a heavy quark Q to the coefficient of the gluon 
scalar-type operator C 9 S as 


C 9 s\q = \ [( a Q + h Q)f+( M 'i m Q’ m i>o) + ( a Q “ h o)f-\ M V m Q^ m i> Q ) 

i=a,b,c 


(53) 


where and (i = a, b, c ) correspond to the contribution of the diagram (i) in Fig. 
They are given as follows: 




A 2 

(— M 2 + m\ + 2-m 2 ) A + — m\ + m|) 

6A 2 ’ 


A 2 


m 2 {(-2M 2 + m\ + 2ml)A ~ 6 m\ml(M 2 + m{ — ml)} 

6?niA 2 


/| 6) (M; mi,m 2 ) = f+\M-,m 2 ,m 1 ) , 
f (p \M;m 1 ,m 2 ) = /i a) (M;m 2 ,m 1 ) , 


with 


A(M; mi,m 2 ) = M 4 — 2 M 2 {m\ + m^) + (m^ — m^) 2 , 




n , -n-i-^+vgh (A > 0) 

< A <°> 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
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In particular, if m\ <U M, m 2 , the above functions are approximated by 


4“’ = 

4* 1 - 


2 ml - M 2 
6(7713 — M 2 ) 2 ’ 
m 2 


3mi(ml — M 2 ) 

4° - 1 


/i 6) ~ 0 , 
/| C) ^ 


6(ml - M 2 ) ’ 

1 


2(m 2 -M 2 ) ’ 


/i c ' ) ~ 0 . 


(62) 


By using the expressions and the tree-level result in Eq. (32), one readily obtains the 
matching condition (35). 
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